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Abstract. We present a unified three-dimensional model of the convection zone and upper atmo- 
t O ■ sphere of the Sun in spherical geometry. In this model, magnetic fields, generated by a helically 
j /^ ' forced dynamo in the convection zone, emerge without the assistance of magnetic buoyancy. We 
use an isothermal equation of state with gravity and density stratification. Recurrent plasmoid 
ejections, which rise through the outer atmosphere, is observed. In addition, the current helicity 
of the small-scale field is transported outwards and form large structures like magnetic clouds. 
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1. Introduction 

H ' 

Usually, magnetic phenomena in the atmosphere of the Sun, e.g., the formation of 
active regions and sunspots, emergence of magnetic fields, and coronal mass ejections, 
are described in terms of magnetic flux tubes. These flux tubes are supposed to form 
y—i . at the tachocline, rise through the convection zone almost undeformed by the effects of 
magnetic buoyancy and reach into the photosphere by creating bipolar regions such as 
On : sunspots. Above the photosphere twisted magnetic fields are observed to form arch-like 
structure. However, there is no clear evidence for the existence of magnetic flux tubes 
inside the deep convection zone. Direct numerical simulations of large-scale dynamos 
suggest that flux tubes are primarily a feature of the kinematic regime, but tend to be less 
pronounced in the nonlinear stage (Kapyla, et al. 2008). In addition, magnetic buoyancy, 
which is usually believed to be the driver of flux tube emergence, can be compensated 
sufficiently by downward pumping resulting from the stratification of turbulence intensity 
in the solar convection zone (Nordlund et al. 1992; Tobias et al. 1998). 

In an earlier work (Warnecke & Brandenburg 2010) we have suggested an alternative 
approach. This is a two-layered approach, where the simulation of a turbulent large- 
scale dynamo in the convection zone is coupled to a simplified solar atmosphere model. 
Magnetic buoyancy does not play an important role in this model. Such a simplified 
model, in Cartesian coordinates, was able to capture some of the observed qualitative 
features. In this work we extend this model further, by going to a spherical coordinate 
system, in the presence of density stratification and gravity. Here we present preliminary 
results from such a study. We find that helical fields, generated by the large-scale dynamo 
below, emerge above the solar surface. We expect such fields to drive flares and coronal 
mass ejections via Lorentz force. 
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2. The Model 

Similar to Warnecke & Brandenburg (2010), a two-layer system is used. We model the 
convection zone starting at r = 0.7 R® to the solar corona till r = 2 R®, where i? is the 
solar radius, used from here on as our unit length. Our simulation domain is a spherical 
shell extending in the 9 (colatitude) direction from 7r/3 to 27r/3 and in the <fi direction 
from to 0.3. This means, the surface in the photosphere, which will be described by 
this simulation, would 600 x 200 Mm 2 large. A helical random force drives the velocity 
in the lower layer. For our model the momentum equation can be written as followed: 

jyT = w (r)f + Vh + g + JxB/p + F vi8c , (2.1) 

where w (r) = \ (l — er ^w)' a P ro f^ e function the connect the two layers, where w 
is the width of the transition, and F v i sc = p _1 V • (2pvS) is the viscous force, S,-j = 
\{Ui-j+Uj-i) — \8ij^-U is the traceless rate-of-strain tensor, semi-colons denote covariant 
differentiation, g = — ^pr-r the gravitational acceleration, h — c 2 lnp is the specific 
pseudo-enthalpy, c s = const is the isothermal sound speed, and / is a forcing function 
that drives turbulence in the interior. The pseudo-enthalpy term emerges from the fact 
that for an isothermal equation of state the pressure is given by p = c 2 p, so the pressure 
gradient force is given by p^ 1 Vp = c 2 V In p = V/i. The continuity equation be written 
in terms of h 

Equations (|2.ip and (|2.2I) are solved together with the induction equation. In order to 
preserve V • B — 0, we write B — V x A in terms of the vector potential A and solve 
the induction equation in the form 

8 A 

^- = UxB + v W 2 A, (2.3) 
at 

For the density we use an initial distribution, where p ~ 1/r 2 . 

The simulation domain is periodic in the azimuthal direction. For the velocity we use 
the stress-free conditions all other boundaries. For the magnetic field we adopt vertical 
field conditions for the r = 2 boundary and perfect conductor conditions for the r = 
0.7 and both 9 boundaries. Time is measured in non-dimensional units r = tu Iuls k{, 
which is the time normalized to the eddy turnover time of the turbulence. We use the 
Pencil Code0, which uses a sixth order centered finite-difference in space and a third- 
order Runge-Kutta scheme in time. See Mitra, Tavakol, Brandenburg & Moss (2009) for 
extension of the Pencil Code to spherical coordinates. 



3. Results 

The forcing gives rise to an a 2 dynamo in the turbulence zone. After a short phase of 
exponential growth, the magnetic field shows opposite polarities in the two hemispheres 
with oscillations and equatorward migration (Mitra et al. 2010). The maximum magnetic 
field at each hemisphere is about 63% of the equipartition value. This is a typical behavior 
of an efficient large-scale dynamo. The magnetic fields emerge through the surface and 
create field line concentrations, which reconnect, separate and rise to the outer boundary 
of the domain. This dynamical evolution is clearly seen in a sequence of field line images 
in Fig.Q] where the field lines of (B)^ are shown as contours of r sin 9(A < f,) c f > and the color 

f http : / /pencil-code . googlecode . com 
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Figure 1. Time series of formation of a plasmoid ejection in spherical coordinates. Contours 
of r sin6{A^)^ are shown together with a color-scale representation of (B^)^,; dark blue stands 
for negative and red for positive values. The contours of r sin (9(^4^,) ^ correspond to field lines of 
(B)<t> in the r6 plane. The dotted horizontal lines show the location of the surface at r = 1 Rq. 
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Figure 2. Times series of what looks a bit like coronal ejections in spherical coordinates. The 
normalized current helicity {J ■ B) $ / (B 2 ) $ is shown in a color-scale representation for different 
times; dark blue stands for negative and red for positive values. The dotted horizontal lines show 
the location of the surface at r = IRq. 

representation stands for (B^)^. Investigating the current helicity, we find a surprising 
result. In the turbulence zone the (J ■ B) $ / (B 2 ) $ is negative in the northern hemisphere 
and positive in south, i.e., has the same sign as the helicity of the external force. Above 
the surface, for each hemisphere, current helicity with sign opposite to the turbulent layer 
is ejected in large patches (Fig. [2]). These structures correlate with reconnection events of 
strong magnetic fields. In Rust (1994) such phenomena have been described as magnetic 
clouds. In order to demonstrate that plasmoid ejection is a recurrent phenomenon, we plot 
the evolution of the ratio (J-B) ( f,/(B 2 ) ( f > versus r and r in Fig.[3J We further find that the 
typical speed of plasmoid ejection is about 0.13 times the rms velocity of the turbulence 
in the interior region, which corresponds to 0.08 of the Alfven speed. The time interval 
between successive ejections is about 100 r. Note that the ejection of magnetic field and 
associated reconnection events are fairly regular, but the formation structures such as 
magnetic clouds are less regular. Note further that the sign of azimuthally averaged 
current helicity in the outer layer is always opposite to that of the turbulence zone. This is 
demonstrated in Fig. [4] where we plot time-series of azimuthally averaged current helicity 
at r = 1.5 Rq. For the northern hemisphere the current helicity (solid black line) and the 
accumulated mean (solid red line) show positive values and for the southern hemisphere 
(dotted lines) negative values. This suggests that, even though the plasmoids observed 
in our simulations are shedding small-scale current helicity of opposite sign to that of 
the large-scale current helicity inside the turbulence zone, outside the turbulence zone 
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Figure 3. Dependence of {J ■ B$ 9)/{B 2 }<f,e versus time r and radius r in terms of the solar 
radius Rq. The left penal show a thin band in 9 in the northern hemisphere, the right one a 
thin band in 6 in the southern hemisphere, both are averaged over 20°-28° latitude. Dark blue 
stands for negative and red for positive values. The dotted horizontal lines show the location of 
the surface at r — 1 Rq . 
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Figure 4. Dependence of (J ■ B) (B 2 ) $ versus time r for a radius of r = 1.5 Rq at ±28° lati- 
tude in arbitrary units. The solid lines stand for the northern hemisphere and the dotted for the 
southern hemisphere. The red colored lines represent accumulated means for each hemisphere. 

the sign of large-scale current helicity has reversed and is now the same as that of the 
small-scale helicity. This may be explained by the action of turbulent magnetic diffusion 
on the magnetic field; see e.g., Brandenburg, Candelaresi, & Chatterjee (2009). 

In summary, it turns out that twisted magnetic fields generated by the helical dynamo 
beneath a spherical surface are able to produce flux emergence in ways that are reminis- 
cent of that found in the Sun. We find phenomena that can be interpreted as recurrent 
plasmoid ejections, which then lead to magnetic clouds further out. A promising exten- 
sion of this work would be to include a Parker-like wind that turns into a supersonic 
flow at sufficiently large radii. 
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